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Abstract 

o , 

The minisuperspace model of a Bianchi Type I universe with compact spatial 
sections is investigated. The classical solutions are brought onto a form were 
the difference between compact and infinite spatial sections are manifest. One 
of the features of the compact case is that it has a non-trivial moduli space. 
The solution space of the compact Bianchi Type I universe is 10 dimensional 
whereas the Kasner solutions only have a 1 dimensional solution space. We 
also include the classical solutions with dust and a cosmological constant. So- 
lutions to the Wheeler-DeWitt equation are obtained in light of the tunneling 

OO ' 

boundary proposal by Vilenkin. Back reaction effects from a simple scalar 
CO ' field are also investigated. 

o 
o 
o 

"o 
cr 

Quantum cosmology was initiated by DeWitt fly in the late sixties. The canonical quan- 
tization procedure which in the twenties was the key to the theory of Quantum mechanics, 
was again used to derive a "Schrodinger equation" for gravity. The idea was that many 
of the cosmological mysteries which had their cause in the initial epoch of the universe, 
could be explained by Quantum cosmology (QC). Through the work of DeWitt, Misner 0, 



I. INTRODUCTION 



Ryan [||-[5| etc., QC had prosperous days in the next decade. Suddenly the interest came 
to a hold, and the number of articles was drasticly reduced. But the interest for QC was 
again renewed in the eighties pioneered through the work of Hartle and Hawking |J. They 
based their work on the boundary condition for a quantum state of the universe. Soon after, 
Vilenkin pj-|T3"[ and Linde proposed their version of a suitable boundary condition. 



Most of the previous work on QC is based on FRW and deSitter universe models although 
some authors have studied anisotropic models (for instance Q). The spatial sections of the 
FRW and deSitter models correspond to globally homogeneous, isotropic, simply connected 



spaces. In recent years a paper by the mathematician Thurston [15| has drawn many a 
physicist's attention. His main conjecture, which still remains to be proven, claims that 
there are essentially 8 types of 3-manifolds. Concerning cosmology and especially QC the 



most interesting are their compact quotients which show interesting properties [|16||. In 
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the nineties several authors investigated locally homogeneous spacetimes with non-trivial 
topologies. Hawking and Turok [17] studied quantum creation of hyperbolic universes in 
the context of the no boundary proposal. Coule and Martin 



in addition to Costa 

and Fagundes [0,0 studied cosmologies with compact hyperbolic spatial sections. We 
also want to mention the work of Fagundes [^TJ already in the early eighties, in which he 
studies compact cosmologies which in the Thurston classification are quotients of K x H 2 . 
In this paper we shall use a locally homogeneous anisotropic model that in the Thurston 
classification has the covering space E 3 . We will not consider in any detail the covering space 
E 3 itself (which is simply connected) but rather one of its compact quotients, the three torus, 
T 3 = M 3 /Z 3 (which is a multiply connected space with a non-trivial fundamental group). 



The possibility of multiply connected spatial sections of the Bianchi metrics |j22| , |23| has again 
renewed the interest for QC. The space T 3 allows a Bianchi Type I metric according to the 
Bianchi classification. 

Using Misner's notation the Bianchi metrics can be written as: 



ds 2 



-N(t) 2 dt 2 + e 2a ^[e 2 ^% jX l X J 



where (3 = diag(/3+ + \/3/3_,/3 + — y/3/3-, — 2f3 + ) and x 1 — Rik^ k where Rik is an SO (3) 
matrix which can be parameterized by its Euler angles (6,4>,ip) i.e. R = e^ K:i e eKl e^ K:i . The 
matrices (ki, /?2, ^3) are the generators of the Lie group SO(3). The invariant 1-forms uj 1 
obey the corresponding Bianchi Type Lie algebra: duj i = —\C % ^ k uji A u k . 

The Einstein-Hilbert action is given byQ S = J M d^x^J —g(R — 2A) + 2 J dM d 3 x\/hK. In 
the case of a spatial homogeneous spacetime like the Bianchi Types, an invariant basis of 
forms on the spatial hypersurfaces can be found. Thus, three of the dimensions in the action 
may be integrated. We assume that we have a topology of a 3-torus, T 3 , and for simplicity's 
sake the volume will be set equal to 1. Choosing compact spatial sections are done for two 
reasons: Firstly because the spatial integration of the action is then finite, but also because 
the solution space of non-isomorphic metrics has a dimension equal to the dimension of the 



(true) phase space p4| 



II. THE MODULI SPACE OF THE TOROIDAL BIANCHI TYPE I UNIVERSE 
AND THE PRINCIPLE OF SYMMETRIC CRITIC ALITY 

Let us briefly review how we construct the torus T 3 as a quotient space of E 3 . It 
is important that we differenciate between the symmetry group that each hypersurface 
possesses and the symmetry group the hypersurfaces possess as spatial sections in the four 



dimensional manifold. For a more thorough investigation on this construction consult |f23 [. 

The full symmetry group of E 3 is 70(3), translations and 0(3) rotations. The Bianchi 
metrics admit a 3 dimensional transitive group of isometries. The 3 dimensional group 
corresponding to the Bianchi Type I Lie algebra is translations in three dimensions, thus 
isomorphic to M 3 . M 3 is a normal Lie subgroup of 10(3) acting simply transitive on the 
spaces M 3 . We will therefore consider the manifold 1R 3 with a simply transitively symmetry 
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group M. 3 defined in the obvious way. We will call this construction M 3 . These spaces have 
exactly the symmetry allowed by the Bianchi Type I Lie algebra. To construct a three torus 
T 3 quotient of M. we can do the following: We find a freely and properly discontinuous 
subgroup T ofQ Sym(IR 3 ), and identify points in M 3 under the action of T. This will ensure 
that the resulting quotient space is a smooth manifold p5| . We are interested in subgroups 
T so that the resulting quotient has the topology of a three torus. These subgroups T can 
be characterized by three vectors a, b, c G IR 3 so that [a, b, c] G GL + (3). These vectors will 
be the generators of the fundamental group of the torus T 3 and points in R 3 are identified 
under the action (which in this case is addition) of the vectors a, b, c. T will then be 
T = {(na + mh + kc) G R 3 |n, m, k G Z}. The moduli space can therefore be described 
by a matrix [a, b, c] G GL + (3) apart from the freedom of discrete modular transformations 
represented by conjugation with respect to matrices in GL + (3,Z). We will not fix the 
gauge with respect to these transformations in this paper. We will simply ignore them. 
The actual parameter space is therefore Teichmuller space. Through the homeomorphism 
GL + (3) = IR 6 x SO (3) we will parametrize the SO(3) sector by the Euler angles (9,phi,ip)- 

In this paper we will follow a slightly different but equivalent point of view. We consider 
a cube C in M 3 . We parametrize the cube by C = {(x, y, z)\0 < x,y, z < 1} and identify 
points on the boundary of the cube so that we have a toroidal topology. The moduli space 
will be determined by the linear mapping /l( x ) — Lx. In the above description the matrix L 
is simply L = [a, b, c] , L G GL + (3). For the Bianchi Type I universe the invariant 1-forms 
can be written locally as u l = K % jdx^ where if* ■ is a constant matrix. After identification 
the coordinates x l become angular-like variables with periodicity 1. 

For a minisuperspace model one always has to ensure that the equations of motion 
obtained from a locally variational principle are equivalent to the Einstein field equations. 
If this is true we say that the Principle of symmetric criticality holds. It has been known 
for a while that the Bianchi Class B models fail to satisfy this principle [p6[| . To prove 
consistency of the minisuperspace model used in this paper we will use a theorem by Torre 
and collaborators p7| , p8i . This theorem states 



The principle of symmetric criticality (PSC) is valid for for any metric field theory 
derivable from a local Lagrangian density if and only if the following 2 conditions 
are satisfied at each point x in the region of spacetime under consideration. 

(1) H«(G,I x )^0 

(2) y/n(KOo = o 

We will explain the symbols as we verify the principle of symmetric criticality in our case. 

G is the symmetry group which is used for reduction. In our case we have G = M. 3 xD 2 
23| . The orbits of G in the spacetime have dimension q, i.e. q = 3. I x is the isotropy group 



for a point x. For a Bianchi Type I universe with toroidal spatial sections each point x has 



2 The symmetry group of M 3 will actually be M 3 xL>2 where D 2 is the dihedral group and x is the 
semi-direct product. D 2 has 4 elements: The unit element and rotations around the 3 axis by an 
angle of it. 
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I x = Z 3 xD 2 c G. H q (G, I x ) is the Lie algebra cohomology of G relative to /. In our case it 
is easy to see that since H 3 (T 3 ) ^ where H 3 (T 3 ) is the deRahm cohomology class of the 
T 3 , we will have H q (G,I x ) 7^ 0. The first condition is therefore fulfilled. V x is the vector 
space of symmetric rank 2 tensors at a point x in the spacetime, and is the vector space 
of /^-invariant symmetric rank 2 tensors at x. (V^)o is the annihilator of V% , i.e. linear 
functions on V which vanish on V 1 . By for instance direct calculation it is not difficult by 
ordinary linear algebra to show that also the second condition holds. Alternatively we can 
use the fact that the group action acts transversally p8|p9| on the bundle of metrics. Thus 



the principle of symmetric criticality holds for the Bianchi Type / universe with toroidal 
spatial topology. 

III. THE KASNER UNIVERSE 



Kasner [32~| solved the Einstein field equations for a Bianchi Type I universe with a 
vanishing energy-momentum tensoif] , T^ v = and spatial sections homeomorphic to the 
Euclidean three-space M 3 . The solution which now bears his name can be written as: 



ds 2 = -dt 2 + t 



tl cos( ^dX 2 + tl^-^dY 2 + tl cos ^ + ^dZ 2 



(1) 



where 7 is an angular variable. If the three (complex) roots of the cubic 



z 3 - ^ = (2) 

are denoted by pi, P2 and p$, the real parts of p±, P2 and p% are equal to the three exponents 
inside the square brackets in eq. [[]. Due to eq. || the following will also be true: 

Pi + P2 + Ps = pi + pl + pI = (3) 

The Kasner solutions are therefore a 1-parameter family of solutions parameterized by 
the angular variable 7. This solution space is therefore usually called the Kasner circle. 

Let us denote the 2-parameter family parameterization of the flat Euclidean 3-space with 
metric 

da 2 (t) = tl cos ^dX 2 + tv^-fodY 2 + tl^+^dZ 2 (4) 

by (1R 3 , da 2 (t)). For a fixed X, Y and Z this family will be volume preserving and homoge- 
neous but it has an anisotropic expansion (in t). 

Despite that Kasner solved the field equations for a vanishing source tensor, the solutions 
for a mixture of dust and a vacuum energy can be brought onto a similar form to that of 
Kasner eq. 



3 Due to the fact that a cosmological constant may be viewed as a vacuum energy, we will treat a 



cosmological constant as a part of the energy- momentum tensor 33]. 
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IV. THE GENERAL SOLUTION OF THE BIANCHI TYPE I UNIVERSE WITH 
T 3 SPATIAL SECTIONS AND WITH DUST AND A COSMOLOGICAL 

CONSTANT 



From 2]| we get the Hamiltonian 



where M is a constant defined by J d 3 x^/hpd us t — 2M. The canonical 1-form is given by: 

= p a da + p+dj3 + + p-d(3- + podQ + p^dcf) + p^dif) 

We have here treated the dust field non-dynamically, i.e. solved the classical equations 
for co-moving dust and thereafter inserted the solutions into the action. The conjugate 
momentum for the dust field is therefore absent in eq. |5] (compare with f3"0,3T|) 



The Einstein field equations are now equivalent to the set of Hamiltonian equations 

-Pi = 7^—, Qi = jr~ ( 6 ) 

dqi dpi 

together with the constraint equation 

H = 

First we will choose the gauge N — 1. From the Hamiltonian equations of motion the fol- 
lowing variables will be constants: The Euler angles (6, <fr, ip) and their conjugated momenta 
{pe,P(j>,Pti}) and the conjugated momenta p±. Let us set 

p + = 4a, p_ = 46. 

and introduce the anisotropy parameter A 2 = a 2 + b 2 . In addition two more constants C± ap- 
pear when the equations for (3± are integrated. The constants (9, 0, ip^pe^p^^p^, A, C+, C_) 
have a clear geometrical meaning. These constants (for A^O) specify a matrixQL e GL + (3) 
such that the regular solid cube is mapped into the spaces (M 3 , da 2 {t)) by the linear map- 
ping Zl( x ) = Lx. The cube is mapped onto a parallelepiped with volume A = det(L). 
Introducing the volume element v = e 3a , the constraint equation yields the equation for v: 

v 2 = 3Av 2 + 3Mv + A 2 (7) 

The constant M will scale as A, so it is more useful to introduce a new constant m by 
M = Am. The volume element v will also scale as A so we introduce the function V(t) by 
v {t) = AV(t). 



4 For an explicit expression of this matrix consult [23] 
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A. The case A > 



A singularity seems unavoidable since for v = 0, v = ±A, the volume element starts 
off at the initial singularity or collapses at the final singularity at a speed A. Our interest 
will be mainly on the expanding solutions. Defining k — ^- and u 2 = 3A we can write the 
solutions as: 

V(t) = — [sinhu;t + KUj(coshujt — 1)] 

where we have chosen V(0) = 0. 

The remaining equations can be integrated in a straight forward manner, and introducing 

, , -(coshcut — 1) 

E + (t) 



sinhajt + klu (cosh, cut — 1) 
the equations for (3± can be solved to yield 

/3 + =isin( 7 -^)lnS + (t) + C + (8) 

0_=±cos( 7 -|)lnE + (t) + CL (9) 

The line-element for a dust-filled Bianchi Type / with a positive cosmological constant is 
therefore 

ds 2 = -dt 2 + \z + {t)l cos ^dX 2 + E+^s^T-i^dy 2 + E+^l^^+i^dZ 2 ] (10) 
whereQ ciX = Ldx and < x, y, z < 1 are "angular" variables. 



B. The case A < 

Let us now consider the case where A is negative. Perhaps the simplest way to obtain 
the solution for A < is to redefine u 2 = 3|A|, and k = ^f, and perform the following 
mapping: 

u i — > iuj 

K I ► — K 

to obtain for the volume function 

V(t) = — [sin cut + kuCI — cos cut)]. 

UJ 



5 Here the notation is a bit misleading(which it often is concerning angular variables). The "exact" 
differentials dX and dx are only exact in the covering space R 3 . On the torus dx can be a globally 
defined 1-form which will be closed but not exact i.e. it corresponds to a non-trivial element in the 
deRahm cohomology class of the torus. 
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We can then define 



E_(f) 



^1 — cos cut) 



smut + ku(1 — COS Ult) 
The resulting line-element is then 

ds 2 = -dt 2 + V(t)* [s_(t)t cos(7) dX 2 + E.{t)l cos ^-^dY 2 + S_(t)i cos(7+ i^rfZ 2 



11^ 



This line element describes a universe that expands at first. After a while the cosmo logical 
term becomes dominant and turn the expanding phase into a contracting one, and it ends 

in a final singularity after an elapsed time t F = - it — arctan 



2 . /W 
3 



C. The case A = 0, and a general form of the line element 

The equations for A = can be solved in a similar manner, and we will therefore just 
write down the result. Defining 



£ (t) 

the line element for A = takes the form 

2 

3 



1 + \mt 



ds 2 



-dt 2 + ( t + -rat 2 
4 



(12) 



Up to a simple rescaling these three line elements are those derived by Saunders ||34||. The 



A = case was also derived by Stephani [ 35] and the m = case by Gr0n [|3(| . Already now 
we can see the resemblance of these three cases and the Kasner line element. However, the 
line element can be brought into a more Kasner-like form by choosing another function N. 
We define a new time variable, r, by: 



dz 



(1 



jmz) 2 — |Az 2 



E: x (r) ,A<0 
EoV) ,A = 

:iHt) ,a>o 



(13) 



where _1 means the inverse function. Inserting this new time variable into the previous line 
elements, we see that all the three cases can be brought into a similar form. Thus we can 
write all the three line elements as: 



ds 2 



dr 2 



((1 - |mr) 2 - |Ar 2 )' 



+ 



(d 



|mr) 2 



At 2 ) 



2 
T3 



r |co S ( 7)rfx 2 + r |cos( 7 -f) dy2 + r |cos( 7+ ^ )ciz 2- 



(14) 
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As we see, the line element is now manifestly cast onto a Kasner form. The solution space 
is now explicitly decomposed into the following sequence: 

SO (3) x R 6 S 1 T w 

1 1 ( 15 ) 

C 1 ± ► R 3 -^l (R 3 ,da 2 (t)) (Rx £,ds 2 ) 

These mappings can be given the following interpretations: 

/l : Deformation of the cube. If we represent the unit cube C by the embedding C = 
{(x'^i = 1, 2, 3|0 < (x i - zj) < 1} for any (x*), then Z L (x) = Lx for L e GL+(3) ^ 
50(3) x M 6 . (See figure[l|) 

s 7 (t) : An explicit 2-parameter parameterization of the Euclidean 3-space with an anisotropic 
metric. Given (X, Y, Z) G R 3 the metric is given by eq. f|. 

K : The Kasner foliation, given explicitly by: K(R 3 , da 2 (t)) = (RxE, ds 2 ) where E = R 3 

and ds 2 = —N(t) 2 dt 2 + V(t)^da 2 (t). The N and V are determined by the Einstein field 
equations and as shown they are given by N(t) = j- — 3 — \ 2 3 v and V(t) = t ■ N(t). 

The solution space is now parameterized by these mappings, and the effect of choosing 
compact spatial sections is now explicitly shown. We have also factorized the solution space 
into a topological (or modular) sector (the 5*0(3) x R 6 sector) and a Kasner sector (the 
Kasner circle). Contributions from the dust and the cosmological terms are also factored 
out. Different matter configurations only affect the last mapping K, all the other degrees of 
freedom affect only the mappings Zl and s T 

Interestingly, the Kasner limit, where A — ► 00 in such a way that lim^ >QO /l(O) = R 3 , 

is well defined. The topological sector becomes degenerate and effectively we have only a S 1 
degree of freedom. On the other hand the FRW limit A — > is apparently badly defined. 
Again the topological sector becomes degenerate in addition to the S^-sector, but we do not 
seem to uniquely reproduce the FRW solutions. 

There is also one more interesting consequence of a non-trivial topology. For the Kasner 
universe (T^ = 0) the special case 7 = 0: 

ds 2 = -dt 2 + t 2 dX 2 + dY 2 + dZ 3 

is through the transformation T = t cosh X, X — t sinh X seen to be equal to the inside of 
the future lightcone in the (T, X) Minkowski space. Thus it can be naturally expanded to a 
flat manifold containing the apparent singularity t = 0. The compact case on the contrary 
can not. Assuming L = As ■ (id) the differential structure of the (T,X) space is that of 
a cone. The singularity t = is now the locus of the cone. It therefore represents a true 
singularity in the C°° structure of the maximal extended space (i.e. the cone including the 
locus). Interestingly we will actually know the true nature of the singularity. Taking the 
limit A — > 00 to obtain the Kasner universe, topological considerations suggest that the 
lightcone in the "flat" (T, X) ought to be identified as one point. The geometry of the space 
7 = is no longer globally flat, it is more like that of a cone. This explains perhaps why 
the "flat space" case 7 = does not evolve as a flat space through the factorization |15[ 



S 



V. THE WHEELER-DEWITT EQUATION 



In the previous sections we have solved the classical equations for a Bianchi Type I 
universe. We will now quantizise the dynamical degrees of freedom according to the Wheeler- 
DeWitt procedure. The constraint equation will then turn into an operator equation on the 
wave function. The resulting equation is the so-called Wheeler-DeWitt (WD) equation. Let 
us also include a real scalar field which later will yield interesting back-reaction effects. 

The WD equation for a Bianchi Type / with a scalar field reads: 



-(-V 2 + V^ = 



where 



-3a 



£a 







V, 



E 



12 V 9a 
--V($)e 3a + AM 



d_ 

da 



+ Vi + 12 



d 2 \ 
<9$ 2 ) 



Here £ represents some of the factor ordering ambiguity and V| 



d 2 



+ 



d 2 

Laplace operator in /5-space. To simplify the expressions we introduce the volume element 
v = e 3a and do the rescaling: f3± = 3/3±, $ = n = |M and V($) = §^($). We 

now drop the tildes to avoid unnecessary writing. The WD equation then turns into (with 



is the usual 



c = i - !) 



dv< 



, 9 
ov 



V 2 



+ fxv + V($)v 2 



(16) 



Since the WD equation does not involve any of the Euler angles, the 5*0(3) sector can be 
treated separately. 



A. The SO(3)-sector 

Since 50(3) = P 3 , this section of superspace is naturally equipped with an elliptic Rie- 
mannian structure. Due to the compactness of 50(3) as a topological space there exists a 
complete countable set of orthogonal functions on this space. Thus the conjugated momenta 
are quantized. We may choose the eigenfunctions in such a way that they are character- 
ized by three quantum numbers: (l,m,m'). These numbers are all integers and obeys: 
—I < m',m < I and I = 0,1,2... The eigenfunctions are not surprisingly the irreducible 
representation matrices of the group 50(3), often written as D® m ,(<j), 8, ip). These functions 
are given by: 

D'LrttM) = e- m<p e-^d^ mf (9) (17) 



where 

A 1 ) (a) _ V l -1 ^ IV! Z " L,iyL ~ " L,iyi Z " L Z " L /-, o\ 

a mm'{ ) , _ , _ , , _ \, K °) 



v (-!)*[(/ + m )!(Z - m)\(l + m')\(l - m')!]* 
n ' { ' V A!(/ + m- 

A 

2l-2X-m'+m / Q \ 2\+m'-m 



X\(l + m - A)!(Z - rri - A)!(A + rri - m)\ 

A 



Q\ At — AA—ift -Tin / Q s 

x(cos-j ( sin-] (19) 
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where the sum is only over integer A which makes sense. The algebraic properties of these 
functions now follow from the group theoretical properties of the Lie group 50(3). We could 
for instance create a "Trace wave function". Every point in P 3 (or element in SO (3)) may 
be represented by a unit rotation vector n in IR 3 and a rotation angle a. Then the "Trace 
wave function" may be defined by: 



sin(Z + \)a 



sin ^ 



These functions have their maximum at the unit element, will be zero where a = Jjrrj 
where n is an integer 1 < n < 21. From a classical quantum description if the Bianchi 
Type I universe was in such a state, the fundamental cube would look as if it were "wob- 
bling" around its identity mapping. It might suddenly do a quantum leap to a different 
orientation. Investigating the "Probability" amplitude l^rl 2 we see that the classical de- 
scription describes a shell model (fig. ||). There are 21 regions where the Bianchi universe 
could be. These regions are separated by forbidden regions. The most probable shell is 
the one containing the identity element. Let us therefore call this the ground state. This 
situation becomes even more evident in the limit I — > oo. It might appear that these wave 
functions allow a geometry change. If a Bianchi Type I universe was born in an "excited" 
state which is meta-stable, it could suddenly cross the barrier to the ground state. The 
state in the neighborhood of the identity element would be the most probable. Due to the 
representation decomposition 

h+h 

i=\h-h\ 



the following will also be true: 



h+h 

^Tr ~ ^Tr ' ^ Tr ~ / ^ * TV 

l=\h-h\ 



Thus the solution space spanned by these Trace wave functions could be generated by the 
two generators: \Pjy and ^jy- These Trace wave functions are of course only a particular set 
of functions which is by no means complete (although orthogonal) but drawn to attention 
only because of the algebraic beauty they possess. The famous orthogonality theorem for the 
matrix elements of irreducible representations will ensure that our functions D^ m ,(4>,6,tp) 
are orthogonal over the Riemannian space P 3 . All of these identities follow from the Lie 
group properties of SO (3) and there are of course a lot more than the ones mentioned here. 
Thus a lot of the properties of the Bianchi Type I are directly related to the structure of 
the superspace, at least in this description. The fact that the superspace has this Lie group 
structure provides a lot of information. 



B. Exact solutions for a zero-mass scalar field 

For a Klein-Gordon field the scalar potential can be written as V($) = A + m 2 $ 2 , where 
A correspond to a cosmo logical constant. For this type of potential exact solutions for the 
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WD equation is very difficult if not impossible to obtain. However, if we are dealing with a 
zero-mass scalar field m = 0, exact solutions are possible to obtain. Firstly we can separate 
out the f3± and $ dependent parts by introducing ^(v, /3±, <3?) = p(p^ e l ( k +P++ k -P-) e m ® . As 
we see k 2 = k 2 h + k 2 _ corresponds to the classical anisotropy parameter A 2 i.e. the volume of 
the fundamental domain. Thus, these particular solutions are planar wave solutions where 
F{v) satisfy the equation: 

v2 ^ + ^ v( ^ + ( k2 + n2 + L iV + Xv2 ) F = ( 20 ) 
This is Whittaker's equation and its solutions can be written as 

F(v) = v a -^W^ L>p (±2i\^v) (21) 

where a — |, L = , p 2 = a 2 — (k 2 + n 2 ) and Wl, p {z) is the Whittaker function. These 
function have an essential singularity at infinity and have a branch point at the origin. 
This is due to the fact that the classical solutions have a physical singularity at vanishing 
3-volume. In the absence of dust /i = 0, the equation |2(J turn into a Bessel equation. Thus 
the solution for \x = can be written as 

F(v) = v a Z p (X^v) 

where Z p (z) is one of the Bessel functions or a linear combination of them. This can also 
be seen from the identity 



W 0>p {2z) = 




The behavior of these solutions are not drasticly altered with the inclusion of a non-zero 
dust term. For simplicity's sake we will set \x = in the rest of this paper. In the further 
investigations the dust term is not essential and because the properties of the Bessel functions 
is better known we will simply drop the dust term. The solutions are therefore (a linear 
combination of) the Bessel functions. 

In the case of a classical state that is not allowed, for instance if we have a negative 
cosmological constant and a large enough v, we would expect an exponential decrease of 
the wave function. For negative A we can write the solutions as a linear combination of 

the modified Bessel functions K p (\X\^v) and I p (\X\^v). However, J p (|A|^t>) ~ e^ v for 
\\\*v 3> \p\, so we have to abandon these solutions if we demand an exponential decrease 
for \X\2v ^> \p\. The actual linear combination is a matter of boundary condition, which we 
will leave for another section. 

C. A non-zero mass scalar field: Harmonic oscillator expansion 

For a non-zero mass scalar field the WD equation is difficult if not impossible to solve 
exactly. 
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Let us try to find solutions of the form: 

oo 



(22) 



n=0 



where ip n {v, 3>) are the eigenfunctions of the harmonic oscillator equation: 



1 d 2 1 2 2 a 
•-— — ■ H — m tr<P 

2 <9$ 2 2 



(23) 



Thus, we can consider the wave equation $ as a path in the Hilbert space which is spanned by 
the basis {\n)} where |n) is the usual normalized eigenstate of the one-dimensional harmonic 
oscillator equation. In a coordinate representation these are given by: 



/mv 



H n (^y/mu) 



where if„ is the ra'th Hermite polynomial. Calculating the derivatives with respect to v: 

"1 1 



(24) 



9 I \ 



-mv $ z 



|n) + ^/ — V W$ 1 77. — 1) 



<9t> 2 



-— — -mt>$ 2 + -m 2 v 2 $ 4 
16 4 4 



mt/) 5 $ 3 |n — 1) 
+ ~y/n(n — l)mf$ 2 |n — 2) 
where we have used the identity: 

H' n = 2nH n - 1 

Using the annihilation and creation operators a and a) defined by: 

1 d 



(25) 



(26) 



a =^ 


mv 

/-* 




jmv 

It* 



V2mv <9$ 
1 d 



(27) 



we can write any power of y/mvQ as 



= (a + a f ) 



Inserting eq. ^2] into the original WD equation for a non-zero mass scalar field, and using 
eq. |2B], ^B] and ^7] we can get a differential equation which the coefficients c n (v) shall 
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satisfy. If we instead introduce a new family of functions defined by d n (v) = v^c n (v) the 
resulting equation may be written as: 

=v 2 d" n + (\v 2 + mv(2n + 1 ) + A;2 + ^- ^- \ n ( n + l )) d n 





1 




-V 




2 






+ 


\( 






1 


+ 
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n(n - l)d' n _ 2 - v^nH- l)(n + 2)< 



n+2 



(28) 



n(n - l)(n - 2)(n - 3)rf„_ 4 + + l)(n + 2)(n + 3)(n + 4)d n+4 

The infinite set of equations divides into two, one odd set of equations and one even set 
of equations, which involve respectively odd or even coefficients only. The Hilbert space is 
therefore divided into two: odd and even states. Now it is time for some approximations. Let 
us investigate the properties of the solutions in the limit v — ► 0. In the further calculations 
we will assume that the wave function ^f(v, $) is independent of <£> in the limit v — ► in 
the following sense: 



Given a compact interval I C R and a 5 > 0. Then there exists a e > so that 
for0<t;< e , \§(v,%)\ <5V$ el. 

This is exactly the right boundary condition for demanding that c n (v) — ► a n (v) as v — > 
where a n (v) is the coefficients of a function independent of <3> (constant in $) in the harmonic 

oscillator expansion. Using the test function / = (^) 4 we see that: 



n=0 

Based on these results we can split the total wave function into an even and an odd part 
\& = ^odd + ^even- We extract the solution for the constant function / by defining a new 

family of coefficients by d 2 n = ^l^' ^n- The equation for £ n will then be 

=v 2 ££ + (\v 2 + mv(4n + 1 ) + k2 + \-^- \< n + \)) £ n 



(29) 



- v (n£n-i ~(n+ \)£'n+?j 

+ ^ (n£n-\ ~{n+ ^)£n+l 

If 13 
+ -ln(n- l)£ n -2 + (n + -)(n+ -)S n+2 

This equation can without difficulty be solved to first order in v. The result is: 

S n = v l \l-\ 1 (An + l)mv + Q(v 2 ) | (30) 
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where I = ~ ± y |g — /c 2 . Assuming £ < 6 the real value of the first order term will be 
negative, i.e. to first order the higher exited modes in the expansion will decay compared 
to the ground state. 



D. Backreaction Effects: The massless case 

The preceding results describe an even wave function in $. The expectation value of the 
scalar field vanishes, but the expectation value of <3> 2 diverges as v — > 0. This will have 
consequences for the effective Hamiltonian in the theory. Let us try to replace the energy 
momentum tensor, T^, by its expectation value < T^ v > using the obtained results. This 
effectively means that the scalar field operator is replaced with its expectation value: 

d 2 

+ m 2 v 2 & 2 



In the massless case the expectation value of the energy momentum tensor only contributes 
with a positive constant n 2 . Keeping the volume of the universe fixed (i.e. keeping A 2 fixed), 
the term n 2 will reduce the radius of the Kasner circle (the "radius" in eq. |2|). Note also 
that the equations ^| will still be true. Since the center of the Kasner circle represents an 
isotropic universe, we see that the inclusion of a massless scalar field effectively reduces the 
anisotropy of the universe. 



E. Backreaction Effects: The Harmonic oscillator expansion 



Using the harmonic oscillator expansion we can try to understand the effect of a non- 
zero mass scalar field on the classical equations. In the harmonic oscillator expansion we 
can write 

oo 



d 2 



i 2 2>Tv2 
+ m V <P 



mv 



^|c„| 2 (2n + l) 



(31) 



n=0 



In our case we shall use only the even part, but unfortuneatly the results obtained can not 
be used as they stand; every sum will diverge. We have to regularize the result so that the 
sums may be performed. Let us write: 

r?(n+|) 



\£J 2 ocv 2Re V 



1 nmv 

V 



where k = Re 



i±v/i- fc2 



and rj is some unknown parameter which represents the uncer- 



tainty in the higher order contribution in the expansion (0). To first order, however, we see 
that this is exactly the previously obtained results. With this replacement the sums may be 
performed and the result is to first order in v: 



even 








^ even) 


\ even 


* even/ 



mv 



277-I 



1 

+ - 

K 



(32) 
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Luckily the unknown parameter r/ does not contribute to lowest order, so the validity of 
the regularization procedure which we used is strengthened. If we assume £ = and we 
associate k 2 with the anisotropy parameter A 2 , we see that we get an effective anisotropy 
parameter A 2 ^ = 2 A 2 + ^. Comparing this to the classical solutions, an increase of the 
effective anisotropy parameter in the Hamiltonian constraint correspond to an increase of 
the speed of the volume element by and a decrease of the anisotropy speed by the same 
amount. 

Even though these considerations were very simple and only approximative, this indicates 
that the effect of a simple scalar field with mass is that it reduces the effective anisotropy. 
Comparing with the zero mass solutions, we see that a zero mass scalar field also has the 
same effect, it increases the anisotropy parameter but reduces the anisotropy. The result in 
the limit v — > for a non-zero mass is thus reasonable. 



VI. BOUNDARY CONDITIONS 

Concerning boundary conditions, quantum cosmology has mostly invoked two different 
boundary conditions. One is that of Hartle-Hawking ||, which is stated in terms of the path 
integral formalism. The Hartle-Hawking wave function ^ hh for a Bianchi Type / universe 
with a positive cosmological constant has been derived by Duncan and Jensen ||37|| . Duncan 



and Jensen also use three-torus as spatial sections, but they do not give any comments 
concerning the SO (3) sector. Apparently the only reason they use compact spatial sections 
is to make the spatial integration in the action finite. In this paper we will however be more 
interested in the tunneling boundary proposal by Vilenkin M. 



A. The Bianchi Type / minisuperspace 

Essential in this discussion is the description of the Bianchi Type / from a minisuperspace 
point of view. If we exclude a scalar field degree of freedom, the minisuperspace will be 
6-dimensional. Topologically we can write the minisuperspace as M 3 x P 3 where P 3 = 
SO (3) and M 3 is the 3-dimensional Minkowski space. As the SO (3)-sector has already been 
discussed, we will only look at the evolution in the Minkowski sector. We will also choose 
the lapse to be unity N = 1. The DeWitt metric on the Minkowski sector will then be: 

ds 2 = e 3a {-da 2 + d(3 2 + + d(3 2 ) (33) 

Doing a conformal transformation we can map this metric onto a Penrose- Carter diagram. 

Mapping the classical solutions (with no dust) into this conformal diagram we get a 
picture like figure |3|. Notice that all the three different cases start off at past light-like 
infinity X~ . They all start off as a Kasner universe (straight line). They end up however, at 
different places. For a positive-valued cosmological constant the universe ends as a deSitter 
universe: moving along the t-axis to the time-like future i + . The Kasner universe (A = 0) 
ends at the future light-like infinity X + , while the case A < the universe collapses again in 
past light-like infinity X~ . The whole boundary except past time-like infinity i~ represents 
singular space-time geometries, and thus according to the tunneling boundary condition, 
is not allowed. However, from the calculations done in the previous chapters, there may 
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be indications that the starting point X~ is somewhat unstable concerning inclusion of a 
scalar field. Even a rather modestly behaved scalar field effectively reduces the speed of the 
universe near the initial singularity to below "the speed of light" . The reduction of the speed 
of a Kasner universe is shown in figure (|. This effectively changes the path's end-points, 
so that the universe starts off at i~ . It ought to be emphasized that this is not enough to 
ensure regularity of the universe at the starting point, but this may indicate that there exists 
a homotopy of time-like paths connecting the Kasner universe with the deSitter solution. 
These different solutions will be universes with different speeds in superspace. One might 
say that we have given the Kasner universe a mass, and the homotopy can be defined by 
the action of the 3-dimensional Lorentz group. 

In the light of this discussion we change the formulation of the tunneling proposal to 
include all universes which may be connected to universes which satisfy the regularity con- 
dition, by the action of the Lorentz group in superspace. 



B. The tunneling wave function 

The tunneling wave function will have a current in superspace which points outwards at 
the singular boundaries of superspace. The conserved current 

J = - ^W*) 

which is analogous to the conserved Klein-Gordon current, states that its zero- component J° 
is a conserved entity. The zeroth component of the Klein-Gordon current can be interpreted 
as a charge, which locally can obtain both negative and positive values. Also the zeroth 
component superspace current can be both negative and positive. In light of the tunneling 
proposal, these correspond to expanding and contracting solutions. 

Using the solutions to the WD-equation for a zero-mass scalar field, we get the "space" 
components of the superspace current to be: 

J ± = - % -v-^(^*d±-$ - #d±#*) = v~f k±\^\ 2 (34) 

At the singular boundary v — + oo, we use the solutions for a positive cosmological constant 
in the big geometry limit. The two possible solutions are the Hankel functions H^^ 2 \v) m 
■^e ±m . The zero-component then becomes: 

jO^ T |A|Va-i)|vl/| 2 (35) 

(2) 

The outgoing modes are therefore the second Hankel function Hp (t>) and the tunneling 
wave functions are are a linear combinations of the functions 

V(v,P ± ,$) = viH ( p 2 \\^v)e i{k+ ^ ++k - p - ) e inq ' (36) 

The solution space of the WD equation will therefore be spanned by the solutions ( |36"D 
multiplied by the solutions from the 5'0(3)-sector D^ m ,. 
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C. Does the Tunneling boundary proposal predict an isotropic universe? 



As we have constructed the Tunneling wave function we might wonder if the wave func- 
tion predict an isotropic or an anisotropic universe. In the classical case we had trouble 
with reconstructing the FRW universe in the limit A — > 0. The FRW universe appeared 
as if it was disconnected from the Bianchi Type / solutions. This picture is drastically 
improved in the Quantum case. If we hold v constant, that is keeping the volume of the 
universe constant, and let A approach zero from a classical point of view we have to let the 
function V(t) approach infinity. For a positive cosmological constant the limit V — > oo 
is indeed a isotropic FRW limit with finite spatial volume. As v and A are independent 
variables we could interpret the A i — ► a FRW limit. Inserting a scalar field we saw that 
it also reduced the anisotropy. The anisotropy was reduced as the ratio k z k +n -z ■ These two 
mechanisms for reducing the anisotropy affect the classical solutions in two different ways. 
Whereas the reduction of A isotropize the universe through the modular and evolutionary 
degrees of freedom, the inclusion of a scalar field reduces the Kasner circle. Thus a scalar 
field could also work for isotropy reduction for the Kasner universe (with infinite spatial 
sections). A — > for the Kasner universe is however more subtle because apparently we 
do not have any clear geometric meaning of the entity A in that case. 

The k 2 enters the wave functions through the order of the Hankel function. The order 
p is given by p 2 = |g — k 2 — n 2 . The factor ordering parameter is assumed to be small, so 
we assume that p is purely imaginary: p = ir where r > 0. A small r means an isotropic 



universe. From the integral representation of the Hankel functions 38 we can write 



T7T 



H%\z) = / e ~ i2COshi - 4rt rft 

n J-oo 

Apparently the wave functions are suppressed by an exponential factor as r grows. Thus the 
probability amplitude is exponentially damped in r: |\I/| 2 ~ e~ rn . We do have to be a bit 
careful however because any linear combination of these wave functions is a solution. We 
might imagine that it had an exponentially growing prefactor in r. Let us try to estimate a 
reasonable prefactor by "normalizing" the wave functions over superspace. Since the variable 
v will behave as a time variable in superspace we could also say that we divide by the "time 
average" of the particular solution. The "normalizing" integrals J dv \^{v)\ 2 diverge for the 
Hankel functions. We will therefore do the following: we Wick-rotate the volume element, 
v i — > iv so that the Hankel functions Hp are "rotated" to the Bessel functions K p . This 
function will under reasonable conditions be square integrable. We estimate the prefactor 
by demanding that the wave function is normalized after a Wick-rotation of the volume 
element. The Bessel functions are related through the following equations 

K p (z) = liei*H®(iz) (37) 
hJ\z) = h£\z) (38) 

We also introduce the integration measure \/\G\ which reflects the geometry of minisuper- 
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space. The DeWitt metric (eq. |33|) will yield \J\G\ = but if we assume^ y/\G\ 
the resulting integral become surprisingly simple p8|: 



lOO A POO 



4 



2 



dv^/\G\\y\ 2 = —e wr dvv K ir {\*v) (39) 

Jo 

' (40) 



Att 2 1 - e- 2 * r 

Using these wave functions after "normalization" the probability amplitude decays as 



e -7TV 



(l-e" 27 "-) (41) 

These wave functions will peak for r = so they predict an isotropic universe. They will 
exponentially decay for larger values of r. Whether this "normalizing" function is the right 
weight function is so far only speculation, but if these functions are reasonable the wave 
functions predict an isotropic universe. 



VII. CONCLUSION AND SUMMARY 

In this paper we have discussed the Bianchi Type I minisuperspace model with compact 
spatial sections. In the classical case this removed some pathologies. The classical solutions 
were all written down in a way so that the solution space was manifestly factorized into 
three parts: a topological sector, a Kasner sector and a matter sector. The Kasner solutions 
were obtained in the limit A — ► oo where one easily could see that the modular sector 
becomes degenerate. Choosing compact spatial sections also gave a better understanding of 
the topology of the 4-space. For instance the apparently flat case 7 = = T^ u of the Kasner 
universe was shown to have a conical singularity at t — 0. 

In the quantum mechanical case since the (true) phase space has a dimension equal to the 
dimension of the solution space, it is easier to give a precise meaning to the WD equation. 
This paper has also shown the necessity of introducing compact spatial sections when one 
wants to quantize the model, and provides a solution to the WD equation for T 3 spatial 
sections. Also a set of solutions in the £0(3) sector have been given which had particularly 
nice algebraic properties. These solutions have a structure which describes a shell structure 
in the SO (3) sector. It was also indicated how these solutions allow geometry changing 
solutions. 

Inserting various matter configurations in general showed isotropization of the universe. 
The resulting tunneling wave functions appeared as if they peaked at small anisotropy. 
Under reasonable considerations we constructed a wave function which had an exponentially 
decaying probability amplitude for increasing anisotropy. This seems to agree with other 
authors [^,0. 



6 Thus these two integration measures coincide iff £ = |. We would also mention that the integral 
is manageable for any power of v (provided that the integral converges). 
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It is quite clear that universe models with compact spatial sections will yield many inter- 
esting and surprising results. Since these possibilities only recently have come to physicist's 
and cosmologist's attention there is a lot of work still to be done. What we have shown 
in this paper is that even at the classical level, this may reveal interesting properties of 
cosmological models. 
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FIG. 3. The evolution of the Bianchi Type / in conformal superspace 
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